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ABSTRACT 



After almost three years of data taking in sky survey mode, the Fermi-LAT 
has detected 7-ray emission toward the Tycho's Supernova Remnant (SNR). The 
Tycho SNR is among the youngest remnants in the Galaxy, originating from a 
Type la Supernova in AD 1572. The 7-ray integral flux from 400 MeV up to 
100 GeV has been measured to be (3.5±l.l stot ±0.7 s2/st ) xlO -9 cm~ 2 s _1 with a 
photon index of 2.3±0.2 stat ±0.1 S2/s t. 

A simple model consistent with TeV, X-ray and radio data is sufficient to 
explain the observed emission as originating from 7r°-decays as a result of cosmic- 
ray acceleration and interaction with the ambient medium. 

Subject headings: acceleration of particles — supernovae: individual (Tycho) - 
radiation mechanisms: non-thermal 



1. Introduction 

Tycho's SNR (SN 1572 or SNR120. 1+01.4) is classified as a Type la (t hermonuclear 



explosion of a white dwarf) based on observations of the light-echo spectrum ( iKrause et al. 
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20081 ). Its expansion has been observed in the radio (VLA, iReynoso et al.lll999l ) and X-rays 
(Chandra. iKatsuda et al.ll2010[ ). The average expansion parameter is v = 0.47 (radio) to 0.52 
(X-rays) and shows strong azimuthal variations. Along 3/4 of the ri m it is consiste nt with 
v = 4/7, expected from a reverse shock developing into r~ 7 ejecta (jChevalierl Il982l ) . while 
toward the east it is much lower, probably due to recent deceleration in a higher density 
medium. The shock speed is approximately 1400 D^ pc km s _1 except toward the east. 



Tycho's distance is not very well constrained. H I absorption studies (ITian fc Leahy 



20111) favor a distan ce on the near side of the Perseus arm (2.5 to 3 kpc). On the other hand 
Hayato et al.l (120101 ) measured the shocked ejecta velocity to be about 4700 km s -1 for Si, S 
and Ar. They concluded that the distance is 3 to 5 kpc on consideration of the measured 
proper motion. 



The radio flux (IKothes et al 



2006 



]) is 40.4 Jy at 1. 4 GHz with a spectral index of 



0.65. The overall X-ray spectrum (IDecourchelle et al.ll200ll ) is dominated by very strong line 



features of Si, S and Fe arising in the shocked ejecta. However the major fraction (60%) 
of the X-ray con tinuum emission co mes from non-thermal synchrotron rather than thermal 
bremsstrahlung ( Warren et al. 20051), and non- imaging in struments indicate that it extends 
to more than 10 keV JFink et al.l 1 1994! ; IPetre et al.lll999h . Recently the VERITAS ground 
based telescope rep orts a TeV detection with a flux of 0.9% of the steady Crab Nebula 
flAcciari et al.lboilh . 



Several arguments point toward a cosmic-ray-modified shock in Tycho. IWarren et al. 



(120051 ) have argued that the contact discontinuity between the ejecta and the shocked ambient 
medium is too close to the blast wave for a compression factor of 4. The other important hint 
comes from the narrow width of the X-ray synchrotron rims, which is probably due to fast 
cooling of the accelerated electrons behind the b last wave, and requires the magnetic field to 
be amplified by accelerated particles to 200 /iG (jCassam-Chena'i et al.ll2007l ). Recent optical 
observations of a fast shoc k in Tycho hav e even provided the first observational indications 
of a cosmic-ray precursor (ILee et al.ll2010l ). 



In this letter we report the GeV detection of Tycho's SNR with the Fermi-LAT . 



2. Observations and Data Reduction 

The Fermi / ~f-ray Space Telescope, launched in June 2008, has operated almost exclu- 
sively in sky-survey mode since 2008 August 4. The principal Fermi scientific instrument, 
the Large Area Telescope (LAT) has a wide field-of-view of 2.4 sr, which allows it to observe 
the whole sky every ~ 3 hr (2 orbits), a large effective area of ~ 8000 cm 2 (on axis at 1 GeV) 
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and a point-spread function (PSF) narrower than 1.0° (for 68% c ontainment) at 1 GeV. For 
a detailed description of the instrument see lAtwood et al. 



We have selected a data set of 34 months in a square region of interest (Rol) of 20° side 
length, centered on the position of the remnant, selecting events from 400 MeV to 100 GeV 
in the PQ_DIFFU SE clasd3- The choice of selecting events above 400 MeV is motivated by 
the broad PSF at low energy («5? at 100 MeV). Moreover the SNR is only 1.4° from the 
Galactic plane, a region mainly dominated by low-energy photons from the Galactic diffuse 
emission, which makes the low-energy determination particularly difficult and subject to 
systematics. 

In addition, to limit contamination from photons produced by cosmic-ray interactions 



in th e upper atmosphere, a 105° cut on the Earth zenith angle has been applied (lAbdo et al. 
2009h . 



3. Analysis and Results 

The 7-ray spectrum from a point-like source coincident with the Tycho SNR has been 
obtained using the Fermi Science Tool gtlike publicly available from the Fermi Science Sup- 
port Center (FSScfl. gtlike has been used in two ways in binned mode: the first method 
consists in fitting all sources in the Rol assuming a simple power law model in the selected 
energy interval; both the integral flux and the photon index of each source within 5 ° from 
Tycho's SNR are kept free in the fit. Only the pulsar PSRJ0007+7303 in the CTAl SNR has 
been fitted with a power law with an exponentia l cut-off. Moreover, all point-like sources 



included in the 1FGL catalog (lAbdo et al.ll2010ah and located within the selected Rol, have 



been accounted for (Fig.[TJleft). The diffuse emission, including the emission from our Galaxy 
and an isotropic component, is modeled according to the rescaled gll_iem_v02 template op- 
timized for the Instrument Response Function (IRF) P6V11. Figure [1] right shows the map 
of the likelihood Test Statistic for E>lGeV of a region 10x10 ° around Tycho. The 

map has been obtained including all the contributions from the backgrounds (point sources 
and diffuse) except Tycho: a bright spot (TS > 25) is visible in the center of the map, which 
corresponds to the Fermi-LAT detection of Tycho. Another peak, about one degree away 



1 http: / /www-glast. slac.stanford.edu/software/IS / 
A http: / /fermi. gsfc.nasa.gov/ssc/ 

3 The Test Statistic is defined as TS = 2(log(Li)-log(Lo)) with Lo the likelihood of the Null-hypothesis 
model (no point source present) as compared to the likelihood of a competitive model, Li. 
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from Tycho, is visible and has been added to the model. 
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Fig. 1. — left: Fermi-LAT count map of the region of 10x10° for photons above lGeV. The 
crosses indicate the pointlike sources included in the model for the fit. A gaussian smoothing 
of 0.3° has been applied, right: Map of the likelihood TS (i.e., map of the TS values obtained 
for a test source that has been stepped through a grid of positions, with the likelihood re- 
optimized at each position) of the same region for photons above lGeV. All contributions 
from known sources and diffuse emission have been included in the source model and so are 
not apparent here except Tycho, which corresponds to the bright spot (TS>25) in the center 
of the image. 



The position of the point-like GeV source is obtained using the tool gtfindsrc selecting 
events above lGeV. The results of the localization based on a maximum likelihood fit are 
R.A. = 6?45 and Dec = 64? 12 (J2000) with a 68% erro r radius of 03°. S ystematic un- 
certainties are estimated to be about 0.006° as stated in lAbdo et al.l ( ]2010ai ): therefore we 
cannot be more precise in localizing the emission within the SNR. The most likely position 
for the 7-ray emission above 1 GeV is shown in Figure |2] as the 95% confidence contour, 
which is in agreement with the radio and X-ray template. Moreover, a study comparing the 
difference in TS between a point-like source at the above mentioned position and the XMM- 
Newton template yields a small TS-difference of 3.3. Since this difference is not significant, 
given the two degrees of freedom (dof) allowed for the position, this confirms that the lo- 
calization of Tycho's remnant in X-rays is fully compatible with the Fermi-LAT data, while 
no further information on size and morphology can be inferred. Concerning the spectral 
properties, the overall fit from 400 MeV to 100 GeV yields a TS of 33 for Tycho, which cor- 
responds to a detection of about 5a for point-like sources with 4 dof. The measured integral 
flux is (3.5±l.l sta t±0.7 sys t) xl0~ 9 cm _2 s _1 with a photon index of 2.3±0.2 stat ±0.1 S3/st . The 
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additional background source at R.A.,Dec = (3?94, 63?94) has a comparable overall flux of 
(3.3±2.0 S fat) x 10 -9 cm~ 2 s _1 with a softer photon index of 2A±0.2 stat ±0.1 syst . Its significance 
is 5a. 




Fig. 2. — Fermi-LAT TS map zoo med in. The green cont ours are 4.5 keV - 5.8 keV con- 
tinuum band from XMM- Newton (jDecourchelle et al.ll200ll ) and the black line denotes the 
95% confidence area for the Fermi-LAT position. 



The second method consists of dividing the entire energy range in logarithmically spaced 
energy bins from 400 MeV to 100 GeV and in fitting all point sources and the scale factors for 
the diffuse components in each bin. The fitting procedure assumes simple power law models 
for all the sources in the Rol and in each energy bin with a photon index arbitrarily fixed 
at 2. The results of this band-by-band fitting procedure are shown in Figure [3] for Tycho. 
In this figure the data points are drawn with statistical uncertainties only; the systematics 
have been indicated with dashed boxes, calculated making the same fit with different IRFs 
(the Monte Carlo-based P6V3 and the in-flight corrected P6V11) and different models of the 
Galactic and isotropic diffuse emission optimized for each IRFs. 



4. Discussion 
4.1. Association 

A point-like GeV 7-ray source is found in the direction of Tycho's SNR based on more 
than 2.5 years of observation with the Fermi-LAT. The positional coincidence between the 




Fig. 3. — Spectrum obtained by evaluating the flux in separate energy bins. The shadowed 
regions indicate the uncertainties due to systematics in the effective area and in the accuracy 
of the modeling of the Galactic diffuse emission. The solid line reproduces the best fit for a 
power-law model. 



Fermi-LAT 7-ray emission and Tycho's SNR suggests that the 7-ray emission is produced 
by shock- accelerated particles within the remnant. We have also considered the possibility 
to have a 7-ray emission from a pulsar, that is the class of most numerous Galactic objects. 
We have excluded this possibility considering tha t Tycho's SNR is kn own to be the result 
of a type la explosion (IRuiz-Lapuente et al.l 12004 iKrause et al.l 120081 ) and that the closest 
radio pulsar in the line of sight, PSR J0026+6320 (about 0.8° from Tycho), has a dispersion 
measure distance of about 13.6 kpc with a spin-down power of approximately 2xl0 32 erg 
s -1 . Moreover, the 7-ray spectrum measured with the Fermi-LA T does not resemb le the 



spectrum of a typical pulsar usually featuring an exponential cutoff (lAbdo et al.ll2010d ). Also 
the chance probability of finding an AGN (the most populated class of LAT sources) within 
the 95% confidence contour has been investigated. The surface density of an AGN with an 
i ntegral flux equival ent to that measured for Tycho or even greater is less than 8xl0~ 3 /deg 2 
(lAbdo et al.l l2010bl ); therefore the chance probability of finding one within 5' is less than 
10" 4 . 



4.2. General constraints 



We conside r the GeV 7-ray production from the bulk of the SNR (not the special 
eastern region). ICassam-Chena'i et al.l ( 120071 ) have analyzed the blast wave of Tycho in the 
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framework of a model that accounts for the dynamical effect of diffusive shock acceleration 
and is broadly consistent with the location of the ejecta, that requires the shock compression 
at the blast wave to be ~6. That model assumes a self-similar evol ution driven by ejecta 
distributed as r -7 , consi st ent with the expansion m easurements of iReynoso et al.l (119991 ) 



and iKatsuda et al.l (120101 ) . iTruelove fc McKed ( 119991 ) have shown that the blast wave radius 
follows the self-similar solution up to twice the time when the reverse shock reaches the 
central plateau. So that solution is valid for the densities considered below up to a distance 
of 4 kpc or so. 

As discussed in the §Introduction, the main uncertainty is the distance to the source. 
Therefore we need to extract the explicit dependence of the observational constraints on 
distance (expressed in kpc) -D kpc . The explosion energy E 5i (in units of 10 51 erg) and 
ambient density are not known precisely e ither. In order to reproduce t he angular size of 
Tycho at the present time (4.26' in the west), ICassam-Chena'i et al.l ( 120071 ) obtain a relation 
between density and distance, assuming E$i = 1. This can be easily generalized to any 
explosion energy, resulting in 

n H = 388D k -; c £ 5 2 lC m- 3 (1) 



T he absence of a clear detectio n of thermal X-ray emission from the shocked ambient 



gas led ICassam-Chena'i et al.l ( 120071 ) to a constraint on density, namely nn < 0.3 cm 3 at a 



distance of 2.8 kpc. That constraint was based on X-ray brightness and it does not depend 
very much on temperature, so it depends little on E$\. Since X-ray brightness is proportional 
to J n 2 dl along the line of sight, the main dependence is as n 2 H D^ vc . Therefore the constraint 
can be generalized to % < 0.5-Dj7 p ° c 5 cm~ 3 . 

We also account for the fact that the total available energy at the blast wave E^ h is 
not the full explosion energy, as some of it is still locked up in the ejecta. In the self-similar 
phase, it can be estimated as 



E 



amb 
51 



0.24n]i 2 DlLlE, 



'kpc-^51 ( 2 ) 

Extrapolating that formula beyond the time at which the reverse shock enters the central 
plateau is reasonable. The transfer of energy to the shocked ambient gas continues, coming 
then from the shocked ejecta. E^ h /E 5 i would be 0.93 at the end of the validity of the self- 
similar phase defined above, so this approximation connects smoothly to the Sedov phase in 
which Elf 3 = E 51 . 

For either bremsstrahlung or 7r° decay which use the same target gas, the predicted 
7-ray flux will be 

F 7 cx / C R^ 5 T b 3nH^k p 2 c oc 0.72f CR E 51 nfD k ^ 2 (3) 
where we have used Eq. [2j 
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Eqs. [T] and [3] define, as a function of distance, a family of solutions in which no decreases 
as -D kp 3 / 2 and E 5 i increases as D^ c 4 . 



4.3. 7-ray emission 

There are three radiation processes potentially responsible for the GeV 7 rays from 
(or in the direct vicinity of) Tycho's SNR: inverse- Compton (IC) scattering on the cosmic 
microwave background (CMB) by relativistic electrons; nonthermal bremsstrahlung by rel- 
ativistic electrons; and 7r°-decay 7 rays re sulting mainly from inelastic collisions between 



relativistic protons and ambient gas nuclei ( iGaisser et al.lll998l ). 



Two different cases are considered in the following: in the first case, called the "nearby" 
scenario, the total energy output of the supernova will be fixed to a standard value of 10 51 erg, 
which results in a distance of 2.78 kpc assuming the maximally allowed value for the ambient 
density of 0.3 cm -3 . In the second case, we place the remnant at a distance of 3.5 kpc and 
calculate the supernova energy and its ambient density to be 2xl0 51 erg and 0.24 cm -3 
respectively. This constitutes the "far" scenario. 

In both cases, the synchrotron flux is constrained by the radio and X-ray data. These 
data imply that a population of shock-accelerated electrons described by a power-law spec- 
trum with a spectral index of 2.2-2.3 and a cut-off energy of 6-7 TeV is the origin of the 
observed synchrotron emission if we a ssume the downstream magn etic field to be ~ 215/xG 



as inferred by X-ray measurements (jCassam-Chenai et al.l 120071 1 . Since the magnetic field 



strength is the only parameter to determine the strength of the IC flux produced by the 
above mentioned population of electrons on the CMB radiation and IR photon fields, the IC 
flux at ~ 1 GeV is inevitably far below the observed value. 

Bremsstrahlung is the only way to possibly fit the data in a leptonic model. The 
flux of bremsstrahlung 7 rays scales with the electron population and n-^. Since the electron 
population cannot be increased too much, otherwise the IC contribution would quickly exceed 
the TeV data, the ambient density needs to be increased and the downstream magnetic field 
has to be decreased to 65 /iG. In the most favorable case of the nearby scenario, the required 
gas density would be riu = 9.7 cm" 3 (using an effective density twice larger behind the 
shock, as prescribed by a Sedov model). This value exceeds the gas density allowed by X-ray 



measurements by a factor of 30 (jCassam-Chenai et al.ll2007l ). Moreover, given the extremely 



high gas density, the energy in the accelerated hadrons needs to be severely reduced to at 
least 1.5 xlO 48 erg which would correspond to an extremely high electron to proton ratio 
(K ep = 0.1). In this case the SNR has to be in the Sedov phase and thus the supernova 
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energy would be 4.4 x 10 51 erg, which is exceptionally high for a type la supernova. For 
all these reasons the bremsstrahlung and IC channel are very unlikely to account for the 
Fermi-LAT measurement. 

On the other hand, the expected 7-ray spectrum of hadronic origin can be calculated 
on the assumption that efficient proton acceleration is taking place at the forward shock 
in Tycho's SNR, whereby the protons acquire the same power-law spectrum with the same 
spectral index as the aforementioned electron spectrum. Only the cut-off energy is much 
higher as protons do not suffer from synchrotron losses as much as electrons do. Since 
the VERITAS spectral measurements (lAcciari et al.ll201ll ) do not indicate any cut-off, we 
estimate the maximum proton energy by equating the acceleration time (assuming Bohm 
diffusion) with the age of Tycho. For = 215 nG , this results in an energy break for 
protons of E P:ma , 



44 D kpc 2 TeV flParizot et aU 12009 ). 



The intensity of this emission depends on the total energy in the accelerated protons 
(and other ions) as well as on the density of the ambient medium %• Even though the 
density is expected to be around 6 % just behind the shock front, the average density seen 
by the shock-accelerated protons over the full emission zone between the blast wave and the 
ejecta is only around 3 wh in the self-similar model. 

Furthermore, the 7-ray emission has been computed assuming that shock acceleration 
is not very efficient and only 10% of the available energy eventually gets transferred into the 
protons/cosmic-rays {E^ h = 0.77i?5i or O.6I-E51 in the "far" or "nearby" case respectively 
from Eqj2]) 

The relevant parameters for the two different cases used here are summarized in Table 
[TJ As shown in Figure H] this conventional hadronic model can explain very well the whole 
7-ray emission from the GeV to the TeV part of the spectrum in a way consistent with all 
the constraints. 



Case 


D kp C 


n H 




Ep,tot 


K e p 


Ep,maic 




[kpc] 


[cm -3 ] 


[10 51 crg] 


[10 50 crg] 


10" 4 


TeV 


Far 


3.50 


0.24 


2.0 


1.50 


4.5 


540 


Nearby 


2.78 


0.30 


1.0 


0.61 


7.0 


340 



Table 1: Parameters used in the spectral energy modeling shown in Figure [H Spectral indices have 
been fixed to 2.3 for both electrons and protons. 
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Fig. 4. — Broadband SED model of Tycho's SNR for the far scenario. 
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5. Conclusions 

A 5cr detection of GeV 7-ray emission from Tycho's SNR is reported. The flux above 
400 MeV is (3.5±l.l stat ±0.7 S2/st ) x 10~ 9 cm^s" 1 and the photon index 2.3±0.2 stat ±0.1 S2/st 
. The measured Fermi spectrum as well as the available radio, X-ray and TeV data can 
be explained by an accelerated proton population which produces 7-ray photons via 7r° 
production and decay. IC emission and bremsstrahlung can account for only a fraction of 
the observed 7-ray flux. 
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